A note on the rare decay of a Higgs boson into photons and a Z boson 
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We have calculated the width of the rare decay H — > 77Z at one-loop level in the standard 
model for Higgs boson masses in the range 115 GeV < tuh < 160 GeV . For this range of Higgs 
boson masses we find that Z boson is predominantly longitudinally polarized, and the photons 
have the same helicity. A comparison of the decay width Y{H — > 77Z) to those of H — > 77 and 
H — > "fZ shows that, for the Higgs boson mass of ran ~ 135 GeV, the ratios of the decay widths 
are F{H -> 77Z)/r(# -> 77) ~ T(H -> j-yZ)/T(H -> 7Z) ~ 10~ 5 - 10 -6 . 

PACS numbers: 13.15.+g, 14.60.Lm, 14.70.Bh, 95.30.Cq 



1. INTRODUCTION 

Due to the lack of a direct coupling between the Z 
boson and the photon in the standard model, the low- 
est order contribution to the decay H — > 77 Z occurs at 
the one-loop level. The Feynman diagrams for this pro- 
cess are identical to those of the crossed channel scatter- 
ing process 77 — > ZH, which was studied by Gounaris, 
Porfyriadis, and Renard The decay H — > 77, which 
also takes place at the one-loop level and is usually viewed 
as a discovery mode for an intermediate mass Higgs boson 
@ , dominates the decay process H — > 77 Z by several or- 
ders of magnitude. However, by imposing kinematic cuts 
on the photons and Z boson in the decay H — > 77Z, 
we may exclude contributions of the back-to-back pho- 
tons. It is therefore in principle possible to distinguish 
the photons from H — ► 77 decay from those arising from 
the decay H — > 77Z. While admittedly very rare, the 
decay H — > 77^, among one loop decay processes, is 
particularly sensitive to top quark couplings. 

In the next section, we discuss the calculations of the 
decay width, the photon invariant mass decay distribu- 
tion, and the decay energy distribution of the Z . This is 
followed by a summary and conclusions. 



2. DECAY WIDTH CALCULATION 

To facilitate our calculations of the amplitudes for the 
process H — ► JjZ we use a generalized non-linear gauge 
fixing condition 3]. In the gauge fixing terms there are 
several parameters which can be set freely without affect- 
ing the values of the measurable quantities. We set the 
values of these parameters in such a way as to minimize 
the number of Feynman diagrams contributing to the 
H — ► 77Z amplitudes. In addition, we use the 't Hooft- 
Feynman gauge, which reduces the propagators for the 
gauge bosons to a simple form — i^^u /(k 2 — to 2 ), where k 



and to are the momentum and the mass of a gauge boson, 
respectively. The resulting Feynman diagrams are iden- 
tical to those of the crossed channel scattering process 
77 — > ZH , which was investigated in the Ref. [l| for the 
standard model and its minimal supersymmetric exten- 
sion. The generic diagrams are drawn |^ in the Fig.^ 
As noted in the Ref. [lj, due to the charge conjugation 
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FIG. 1: The generic Feynman diagrams for the decay process H — > 
77Z are shown. For the charged fermion / in the Figs.^a) and 
^b), we include only the top quark. However, in the Fig. 0c), the 
Z boson exchange, all charged fermions / are included in the loop. 

properties of the gauge boson couplings there are no W 
bosons in the loops and we need only consider contribu- 
tions from charged fermion loops. In the box diagram, 
Fig.^a), the coupling of the Higgs boson to the fermion 
in the loop is proportional to the fermion mass and we 
include only the top quark contribution. The same thing 
is true in the triangle diagram with the neutral Goldstone 
boson G° exchange, Fig.^b), and, again we include only 
the top quark contribution. 

In the Z boson exchange diagram, Fig.^c), we en- 
counter an anomalous contribution |FJ, whose cancella- 
tion requires the inclusion of all charged fermions of a 
given generation. After the cancellation of the anoma- 
lous contribution, which is independent of the to/, the 
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mass of the charged fermion in the loop, the contribu- 
tion of this diagram to the dominant helicity amplitudes 
is proportional to the m 2 . As mentioned above, the con- 
tributions of Fig.nja) and Fig.^b) to any helicity ampli- 
tude is proportional to the m 2 We, therefore, could 
have included only the charged fermions of the third gen- 
eration in the evaluation of diagram Fig.^c) but chose 
to include all three generations in this diagram. 

For the decay H — > 77Z, the helicity amplitudes 
-4aiA 2 a z , where Ai and A2 are the helicities of the pho- 
tons and Xz is the helicity of the Z boson, are not in- 
dependent. Due to Bose symmetry and CP invariance, 
there are only four independent helicity amplitudes |l(, 

which we choose to be -4++o, -4+++, A-\ , and A^ o- 

Any other helicity amplitude is expressible in terms of 
these four amplitudes. For instance, CP invariance re- 
lates -4--o to the basic set as 

A—o = -A+ +0 . (1) 

Our numerical calculation of the amplitudes indicates 
that, among the four basis amplitudes, the helicity am- 
plitude .4++o is by far the largest. As a consequence, we 
need only consider the amplitudes -4++o and A o- 

Our calculations were performed in a semiautomatic 
method [6|. Specifically, we used the FeynArts pack- 
age 0, to generate the helicity amplitudes in terms of 
the tensor loop integrals. Using the algebraic manipula- 
tion program FORM |8(, these tensor integrals were con- 
verted to the conventions of Passarino-Veltman [9| and 
expressed in terms of scalar integrals. The scalar inte- 
grals were calculated using code written following the 
approach of 't Hooft-Veltman 0] an d Denner-Nierste- 
Scharf [Tlj. As a partial test of our calculations, we com- 
pared the numerical results of our routines for the scalar 
integrals with those of the FORTRAN codes LOOP 0] and 
the ff 0] . While these two packages did not give numer- 
ically stable results for all regions of the parameter space, 
where comparisons were possible we obtained essentially 
identical numerical results for our scalar integrals and 
those of the loop and/or ff. We also checked the gauge 
invariance of the helicity amplitudes for each of the final 
state photons. 

To facilitate the discrimination of H — > 77Z from 
H ~* 77 an d H — > 7Z and account for possible experi- 
mental limitations, we imposed a variety of cuts on the 
following photon and Z boson variables: \py\, \py\, \pz\, 
m*y, m 2 z , m^, z , 7 y, 6 lZ , and 9y Z - Here, p y , py, and 
pz are the 3-momenta of the photons and the Z boson, 
respectively, in the center of mass of the Higgs boson 
and 9 7 y, 0~f Z , and 9y Z are the various angles between 
the 3-momenta, p 7 , py, and p z . The invariant mass 
variables are m 2 ^ = (p 7 + py) 2 , m 2 z = (p 7 + p z ) 2 , 
and m 2 , z = (py +p z ) 2 . Here, p 7 , py, and p z are the 
4-momenta of the photons and the Z boson, respectively. 

For our calculations of the decay width T(H — > 77Z), 
invariant mass distribution dT(H — > r y"fZ)/dm~ f y, and 



the Z boson energy distribution dT(H — > jjZ)/dE z , we 
choose the following set of cuts: 

|P 7 |cut = liVlcut = |pz|cut = |p|cut, (2) 

(m 77 ') cut = (m jZ ) cnt = (my z) cut = m cu t, (3) 

(#77')cut = (&yz)cut = (Qyz)cut = 6 cut • (4) 

These cuts facilitate the experimental tagging of the pho- 
tons and the Z boson. They provide minimum opening 
angles between the photons and the Z boson, exclude 
contributions of the back-to-back photons, and also avoid 
any numerical instability in the calculations. The cuts 
help discriminate the non-back-to-back photon pairs of 
the decay H — > "f"fZ from the back-to-back 77 pairs in 
the decay H — > 77. In principle, all the photons of the 
decays H — > 77, H — » jZ, and H — > 77Z can be identi- 
fied. 

In Fig.|5J we show the result of the calculation for the 
decay width T(H — > 77Z) as a function of the Higgs 
boson mass m#. For purposes of comparison, in this 
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FIG. 2: (Color online) The decay widths as function of mjy for 
three decay modes of the Higgs boson are shown. The solid line is 
T(H — > 77Z), the dashed line is T(H — » ~fZ), and the dotdashed 
line is F(H — > 77). The cuts imposed on F(H — > 77.Z) are |p| C ut = 
5 GeV, m cut = 10 GeV, and cu t = fj • 

figure we also include the decay widths T(H — » 77) 
and T(H — » jZ) Ibi, w hich are calculated using the 
hdecay package [jja, Il6| . As it is clear from this fig- 
ure, for the Higgs boson masses that we are consider- 
ing, the decay width T(H — > 77Z) is several orders of 
magnitude smaller than the decay widths T(H — * 77) 
and T(H — > jZ). For instance, for a Higgs boson mass 
of ran ~ 135 GeV, the ratios of the decay widths are 
T(H -» ~ryZ)/T(H -> 77) ~ T(H -> ^Z)/T(H -> 
7Z) ~ 10~ 6 . There are no cuts on the decay widths 
T(H — > 77) and T(H — > jZ). The cuts we imposed, 
of course, decrease the value of T(H — > jjZ). In ad- 
dition, there is the suppression from three-body phase 
space, and from the higher order in a. However, there are 
other differences in the amplitudes for these decay modes, 
which account for their relative sizes. For instance, in the 
case of H — > 77, the decay amplitude receives contribu- 
tions from charged fermion loops as well as a substantial 
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contribution from W boson loops, whereas in the decay 
H — ► 77.Z, there are no W boson loop contributions and 
the anomalous triangle diagram, Fig.^c), must be in- 
cluded. As a result of these differences, the simple power 
counting method for estimating the size of the ratio of 
the decay widths T(H — > jjZ)/T(H — > 77) is rather 
unreliable. 

In Fig. we show the invariant mass distribution 
dT(H — > r yjZ)/dm- ry i as function of m 77 ' and in Fig.|31 
we show the energy distribution dT(H — ► jjZ)/dEz as 
function of the Z boson energy Ez- 
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FIG. 3: (Color online) The invariant mass distributions F(H — > 
■yyZ)/dm 1 ^i as function m 7 y , the invariant mass of the final pho- 
tons, for Higgs masses of rrijj = 115, 125, 135, 145, and 160 GeV 
are shown. The cuts imposed are the same as those for the total 
width V(H -> 77Z) of the Fig.|2| 
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FIG. 4: (Color online) The energy distributions F(H — > 
■yyZ)/dEz, as function of the Z boson energy Ez, for the Higgs 
masses of Fig.[3]are shown. The cuts imposed are the same as those 
for the total width T(H -> 77Z) of the Fig.|2| 

To further investigate the effect of cuts on the decay 
width T(H -> 77Z), in Fig.EJwe plot T(H -> 77Z) as 
function of m#, for several different choices of |p| C ut, 
TO cu t, and 9 CU f In the Fig.|3(a), all three cut parameters 
\p I cut, w cu t and 9 cut are varied simultaneously. This fig- 
ure captures the combined effect of all three cuts on the 



decay width. The effects of the individual cuts on the de- 
cay width are illustrated Figs.|5J(b)-(d). In Fig.[5{b), the 
invariant mass and angular cuts are held fixed and the 3- 
momentum cuts are allowed to vary, in Fig.[S](c) we hold 
the 3-momentum and the angular cuts fixed and vary the 
invariant mass cut, and in Fig.[5](d) the 3-momentum and 
invariant mass cuts arc held fixed and the angular cut is 
varied. From Figs .[51(b)- (d) it is clear that the cuts have 
more effect on the decay width for the low Higgs bo- 
son mass. Fig. |3(d) illustrates that the angular cuts only 
start to show any noticeable effect on the decay width T 
for 6» cut > 7r/6. 



3. SUMMARY AND CONCLUSIONS 

Compared to the one loop decays H — > 77 and H 
jZ, the decay H — > jjZ is highly suppressed in the 
standard model. In one sense this is unfortunate because 
the physics of this process has some interesting features. 
Among these are the absence of W boson contributions in 
any of the loops. The amplitudes are dominated by top 
quark loops and therefore sensitive to top-Z couplings. 
Additionally, there is the presence of an anomalous vertex 
in the s-channel Z exchange diagram, which might be 
studied were it not for the constraints of Yang's theorem 
|l7j . Because Yang's theorem forbids the physical decay 
Z — » 77, the s-channel pole cancels in the mass region 
we consider, 115 GeV < m# < 160 GeV. This, together 
with the anomaly cancellation, serves to further reduce 
the decay width. 

Perhaps the most striking feature of this process is the 
simplicity of the decay amplitude hclicity structure. Al- 
though not apparent at the outset, the calculation shows 
that there are only two helicity amplitudes of any impor- 
tance, -4++o = — A o- As a consequence, the photons 

in H — > 77Z always have the same helicity and the Z is 
always longitudinally polarized. 

In conclusion, we have found that T(H — > 77Z) is very 
small compared to viable one loop discovery modes of the 
Higgs boson and that the suppression is greater than sim- 
ple phase space and coupling constant accounting might 
suggest. It goes without saying that the detection of the 
decay mode H — > 77Z will be extremely difficult if it 
occurs at the standard model level and that any signal 
in this channel is very likely evidence of new physics. 
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FIG. 5: (Color online) Fig.0(a) shows F{H -*■ 77Z) for the cuts (|p| cu t, m cu t,0cut)=(l GeV, 1 GeV, tt/48) (solid), (5 GeV, 5 
GeV, tt/6) (dashed), (10 GeV, 10 GeV, vr/3) (dotdash). In Fig-El(b), m cut = 1 GeV, 6> cut = tt/48 and [p | cut = 1 GeV (solid), 5 
GeV (dashed), 10 GeV (dotdash); in Fig,'E](c), \p [ cut = 1 GeV, cu t = tt/48, and m cut = 1 GeV (solid), 5 GeV (blue), 10 GeV 
(dotdash); and in Fig.|3(d), \p | cut = 1 GeV, m cu t— 1 GeV, and # cut =7r/48 (solid), 7r/6 (blue ), and tt/3 (dotdash). 
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